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ABSTRACT 

Phase perturbations due to inclined surface magnetic field of active region strength are 
calculated numerically in quiet Sun and simple sunspot models in order to estimate and 
compare the direct and indirect (thermal) effects of the fields on helioseismic waves, 
ft is found that the largest direct effects occur in highly inclined field characteristic 
of penumbrae, and scale roughly linearly with magnetic field strength. The combined 
effects of sunspot magnetic and thermal anomalies typically yield negative travel-time 
perturbations in penumbrae. Travel-time shifts in umbrae depend on details of how 
the thermal and density structure differs from the quiet Sun. The combined shifts 
are generally not well approximated by the sum of the thermal and magnetic effects 
applied separately except at low field strengths of around 1 kG or less, or if the thermal 
shift is small. A useful rule-of-thumb appears to be that travel-time perturbations in 
umbrae are predominantly thermal, whereas in penumbrae they are mostly magnetic. 
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1 INTRODUCTION 



Time-Distance helioseismology (TD; Duvall et al 



Helioseismic Holography (HH; Lindsey & Braun 



1993) and 



2000) are 



prominent examples of local helioseismic techniques which 
rely on accurate identification of oscillation phase in order to 
probe the solar interior. Active regions, especially sunspots, 
have been particular targets of these methods in an ongo- 
ing quest to better understand the Sun's magnetic activ- 



ity (|Kosovichev, Duvall fc Scherre r 2000; Braun & Lindsey 



|2000 1. However, Lindsey & Braun (2005) and Schunker et al 
( 2005 1 draw particular attention to the "showerglass" effect, 



produced by surface phase perturbations, which obscures 
the subsurface holographic image, particularly in penum- 
brae. Couvidat & Rajaguru (20071 also identify "surface 
magnetic contamination" of TD data in penumbrae, and 



Braun & Birch ( 2008 1 emphasize the frequency dependence 



of travel times. 



Recently, Cally ( 2008 1 employed exact solutions for 



waves in a two-dimensional (2D) isothermal model with uni- 
form inclined magnetic field to firmly identify near-surface 
phase shifts in helioseismic waves. These occur in particu- 
lar at the upper turning point of fast (magnetically domi- 
nated) waves, and at the mode conversion level z oq (where 
the Alfven speed a and the sound speed c coincide). This 
is a purely magnetic effect; no thermal perturbations were 
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introduced. Such direct magnetic effects are clearly prime 
candidates for explaining the observed surface anomalies. 

Unfortunately, the exact solution approach is restricted 
to the isothermal case, and to 2D, and so is more a proof of 
concept rather than a directly applicable model in the real- 
istic solar context. In this paper we relax these restrictions, 
and quantitatively calculate the surface phase shifts due to 
both 3D magnetic and ID thermal effects, individually and 
in combination. 



1.1 Effects of Magnetic and Thermal Anomalies 

It is common practice in local helioseismology to linearize 
the effects on the travel time perturbations. For example, 
|Kosovichev fc Duvall| ( |1997[ ) propose the following expres- 
sion for the travel-time perturbation in the ray approxima- 
tion 



At(x!,X 2 ) 

n-U 



1 Sc SiO c U! r V v 1 
1 1 + - 



I'd 



(a-k) 2 

c 2 k 2 



ds , 



(1) 



where T is the ray path joining the two surface points x\ 
and X2, h is the unit vector along V, c and a are the sound 
and Alfven speeds, a is the Alfven velocity, U is the un- 
derlying flow speed, u is the frequency, k is the wavevector 
with k = |k| the wavenumber, u c is the acoustic cut-off fre- 
quency, and v p = uj/k is the phase speed. The magnetic 
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field is being treated as a perturbation to the background 
non-magnetic state (hence the lack of a 5a term), with 5c 
and 5uj c characterizing the thermal perturbation. Fermat's 
Principle - that the travel time varies only at second order 
with perturbations to the path F - is invoked to avoid a 5T 
contribution. Although lacking an accounting of the effect of 
the Brunt- Vaisala (buoyancy) frequency, and incorporating 
a description of the magnetic influence which is valid only 
if a < c, equation illustrates how At is perceived as 
depending additively on flow (first term), thermal (second 
and third terms), and magnetic (final term) effects. 

However, Weinberg ( 1962 1 proves that Fermat's Prin- 



ciple applies within ray theory only if the dispersion func- 
tion D(tj,k) is homogeneous in lu and k, i.e., D(au),a\i) — 
a m _D(aj,k) for some m, which is not the case if acoustic 
cut-off or Brunt- Vaisala frequencies are included (see also 



Barnes & Cally 2001 1. Furthermore, classical ray theory 



breaks down in the neighbourhood of the mode conversion 
layer z « z cq , where rays are split and their subsequent 
paths radically altered ( jSchunker fe Cally |2006l |Cally|2007| ) , 
again invalidating Fermat's Principle. Surface magnetism in 
active regions certainly cannot be regarded as a mere linear 
perturbation. 

In practice, only flow and thermal perturbations have 
been incorporated into TD inversions up till now. Recently 
though, it has become clear that direct magnetic effects, 
through both magnetic perturbations to the wave speeds 
and to the ray paths, are crucial, indeed dominant, in the 
surface layers of sunspots ([Cameron, Gizon, fc Duvall|2008| 



Moradi & Cally 



2008 



Moradi, Hanasoge, & Cally! 2008 



In particular, based on sophisticated numerical simulations, 
[Cameron, Gizon, fc Duvafll conclude that "constraining the 
sunspot model with helioseismology is only possible because 
the direct effect of the magnetic field on the waves has been 
fully taken into account" . 

A further issue with current practice arises from equa- 
tion |TJ or similar: whereas inversion is based on such equa- 
tions, the "observed" At, which comprises the data being 
inverted, derives from measured phase shifts. However, as 
|Cally| ( |2008| ) demonstrates, phase is not continuous along 
a ray path, but instead suffers significant jumps at turning 
and mode conversion points. It may be hoped that, although 
total travel time measured through phase will certainly not 
agree with that found by integrating phase speed along ray 
paths, the perturbations may be tolerably accurate. But this 
is by no means certain when strong surface magnetic fields 
are involved. 

Mean one-way travel-time perturbations in the range 



—50 s < At < 40 s are reported by Couvidat, Birch & Koso- 
|vichev| ( |2006[ ), with negative At, for deep, comparatively 
\ow-l modes, and positive At for shallow high-£ waves. Using 
ridge-filtering rather than phase-speed filtering, |Couvidat fc| 
|Rajaguru| pTJ07| and |Braun fc Birch] (|2008f find that posi- 
tive travel-time perturbations are more prevalent for shallow 



1 Both | Moradi fc Ca lly ( 2008 ) and Robert Cameron (private 
communication) find that travel-time perturbation magnitudes 
are greatly reduced in most cases if direct magnetic effects are ar- 
tificially "turned off" whilst retaining the thermal perturbations 
in sunspot models, strongly suggesting that the direct effects gen- 
erally dominate. 
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Figure 1. The sound speed profiles near the surface in the 
MMMS models with (top to bottom) kG, 1 kG, 2 kG, and 
3 kG magnetic fields respectively. 



waves in sunspot umbrae, with negative shifts more charac- 
teristic of penumbrae. However, they show that such details 
are very sensitive to the phase-speed and frequency filters 
used, with more extensive negative-Ar regions at higher fre- 
quencies (4.5 mHz). 

With these points understood, the specific questions to 
be addressed here are: 

(i) What "travel-time" shifts does the magnetic field 
alone produce (always bearing in mind that these are con- 
ventionally calculated from the observed phase shifts using 
At = — u>~ 1 A<p, and therefore are not strictly travel-time 
shifts at all when there are phase jumps)? 

(ii) What travel-time shifts does the thermal structure of 
a sunspot produce, in the absence of magnetic field? 

(iii) What is At due to the combined thermal and mag- 
netic anomalies, and is it a simple sum of the individual 
effects? 



2 MODEL 



The model is essentially that of Cally & Goossens (2008) 



The well-known quiet Sun Model S of |Christen scn-Dalsgaard 
|et al.| ( |1996[ ) is supplemented by a uniform magnetic field 
B = B (sin 6 cos <j>, sin 8 sin <j>, cos 9) inclined at angle 6 to 
the vertical and angle <f> to the x-z Cartesian plane in which 
the wave is defined to propagate. Since the field is uniform, 
there is no net magnetic force on the medium, so the model 
is formally consistent. The 2D case corresponds to 4> = 0. 
An infinite isothermal atmosphere is placed above Model S 
(z > 0.5 Mm). The thermal anomaly is crudely accounted for 
in the Magnetically Modifie d Model S (MMMS, see figure [TJ 
of |Schunker fc C ally (2006, Appendix A), which reduces the 
Model S gas pressure by the constant magnetic pressure, and 
hence leaves the density unchanged (except in the overlying 
atmosphere where this pressure balance is untenable). The 
MMMS models take field strength B as a parameter, and 
reduce the temperature near the surface compared to Model 
S, with greater cooling as B increases. Although ad hoc, 
the MMMS models mimic the crucial cool surface layer of 
a real sunspot, and will be adequate to test the relative 
roles of the thermal and magnetic anomalies. An alternate 
thermal model which allows for the atmosphere to compress 
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downward with decreasing temperature is briefly discussed 
in Section 13731 

It has become abundantly clear in recent years that 
sunspots are a very shallow phenomenon, with most of their 
helioseismic impact being restricted to within a few hun- 
dred kilometres of the surface. The density and Alfven scale 
heights are similarly short. This is to be compared with hor- 
izontal distances of several megametres or more over which 
the field and thermal structure vary. It is therefore appropri- 
ate to explore horizontally homogeneous models in the first 
instance, especially when concerned with waves which ap- 
proach the surface nearly vertically, and resolve the umbra 
and penumbra. 

The approach is that of a scattering experiment. We en- 
visage an infinite horizontal driving plane at depth z — Zb < 
injecting fast magnetoacoustic waves with exp[i(k x x—ui t)] 
horizontal and temporal dependence into the overlying near- 
photosphere. The horizontal wavenumber k x is chosen so 
that the wave's lower turning point would be at a prescribed 
depth, Z\ < Zh, if the model were not truncated at z\,. Typ- 
ically, Zb is placed deep enough that a <C c, and the fast 
and slow magnetoacoustic waves are well-decoupled. Since 
the sound speed dominates in this region, fast waves are 
overwhelmingly acoustic, with at most only minor magnetic 
characteristics. The upgoing fast (acoustic) wave is partially 
converted to a fast (magnetic) wave near z eq : mode conver- 
sion in the 2D case, where the Alfven wave decouples, is ex- 
tensively explained using generalized ray theory in |Schunker| 



fe Cally (20061 and Cally (20071. Above z cq it refracts from 



the steep Alfven speed gradient and propagates back down- 
ward. The remainder of the incident wave's energy is trans- 
mitted as a nearly field-aligned slow (acoustic) wave, and, in 
the 3D case, as an Alfven wave. On its way back downward, 
the fast (magnetic) wave again meets a — c and is again 
partially converted to a fast (acoustic) wave, partially to an 
Alfven wave, and partially (transmitted) to a slow (mag- 
netic) wave. It is the downgoing fast wave in z < z cq that 
is of prime interest here, as it is this wave that re-enters 
the helioseismic wave field. It skips in the usual manner, 
and may be observed and interrogated at the surface in the 
surrounding quiet Sun. 

There is another route to the skipping fast wave though. 
If the wave frequency u is less than the ramp-reduced acous- 
tic cutoff frequency lu c cos#, the transmitted slow (acoustic) 
wave in z > z cq will reflect, pass through a = c, and par- 
tially transmit to a fast (acoustic) wave in z < z cq , where it 
will join the fast wave which took the first route. 

We are interested in both the downgoing fast wave's 
energy flux (relative to that of the original injected wave) 
TZ € [0, 1] and its phase. Using numerical solution of the full 
linearized 6 th -order wave equations and a least-squares fit 
of the calculated acoustic field $ = p 1//2 c 2 V- £ to the stan- 
dard WKB acoustic solution, * = C+ /+ (z) + C_ f-(z), the 
upgoing and downgoing coefficients C+ and C_ are found 
(see the appendix for details). Here, p is the density and £ 
is the plasma displacement. The fit is performed over a se- 
lected interval Zb < z < z s chosen to sit comfortably within 
the acoustic cavity, and its accuracy is normally excellent 
in practice. (It does become more difficult though, and the 
results slightly more uncertain, for a few of the very shal- 
low waves with lower turning point Z\ = —3 Mm explored 
in Section |3.2.2| but the uncertainties are more matters of 



detail than they are fundamental.) Then, having defined /+ 
and /_ to be equal at Zb, the reflection coefficient and total 
phase change at Zb are 



11 = 



C- 



C*4 



and 



8p — ~ arg 



(2) 



The raw Sip is not very useful in this form. It can be al- 
tered by simply moving Zb up or down. However, the signif- 
icant quantity is the difference between Sip in the magnetic 
and non-magnetic cases, Aip — Sip — Sipo, both calculated 
with the same Z\. This is independent of 2b, provided it is 
well within the c>a regime. Unlike [Caliyl ( |2008t , we do not 
calculate the individual phase jumps at the conversion and 
top turning points, but only the total perturbation through 
the surface layers. The total inferred phase travel time per- 
turbation, as used in TD, is then At — —to" 1 Aip: phase 
advance Atp > corresponds to negative At, i.e., faster 
propagation. 

Three boundary conditions are specified at the top, Zt, 
situated in the overlying isothermal atmosphere: that the 
fast wave is evanescent rather than exponentially growing 
with height; that the slow wave is either upwardly prop- 
agating or evanescent, according as co exceeds cj c cos# or 
vice versa; and that the Alfven wave is outgoing. |Cally fe| 
Goossens ( 2008 ) present full details of how these conditions 



are imposed. In the spirit of a scattering experiment de- 
signed to probe the surface layers, we do not attempt to 
model the chromosphere. The isothermal layer, in which we 
have exact solutions, is a convenient "observation region" 
in which the transmitted waves may be examined. Further 
reflection off the rise from temperature minimum to chro- 
mospheric temperatures, or off the transition region, as well 
as radiative effects, will produce somewhat altered results, 
but they are outside our focus of interest. 

The mathematical problem is completed by the further 
imposition of three boundary conditions at the bottom, Zb: 
fast, slow and Alfven waves are outgoing, and an incoming 
fast (acoustic) wave is imposed. No restriction is placed on 
the resulting outgoing fast wave, of relative flux 1Z. 



3 RESULTS 

3.1 Thermal Anomaly Only 

Figure [2] displays the travel-time perturbations for the 
MMMS models with thermal perturbations corresponding to 
S e ff = 1, 2, and 3 kG effective magnetic fields, but with the 
direct magnetic effects on the waves suppressed. For compar- 
ison with the observed travel-time perturbations ( jCouvidat^] 
Birch & Kosovichev 2006 1 mentioned earlier, note that those 
presented here are two-way perturbations; i.e., the wave has 
travelled both up and down from and to zt, and so is dou- 
ble the one-way shifts. The cooler models have produced 
increased travel times, with greater effect for larger B e s, 
and lesser at low frequencies, at which the wave turns over 
much lower, and hence does not substantially sample the 
surface thermal perturbation. This is a feature of the MMMS 
models, in which temperature is reduced but density is un- 
changed from model S. Hence, the acoustic cutoff frequency 
uj c — (c/2H)yl — 2H' (where H is the density scale height) 
is also reduced, and if anything the wave reaches higher be- 
fore being reflected. The combination of lower sound speed 
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Table 1. Spherical harmonic degree I for the adopted frequencies 
and lower turning depths. 





3 mHz 


4 mHz 


5 mHz 


-3 Mm 


643 


936 


1213 


-5 Mm 


521 


733 


937 


-10 Mm 


341 


465 


586 


-15 Mm 


274 


370 


465 



3.2.1 Vertical Magnetic Field 

Near-vertical magnetic field of around 3 kG or more may 
be found in large sunspot umbrae. Figure [3] presents the in- 
ferred travel-time shift At = — u~ Acp as a function of field 
strength for several frequencies and lower turning depths. 
(Note though that At is simply a measure of phase shift Acp, 
and is not necessarily a true travel-time shift.) As expected, 
At diminishes with reducing B, though the limit cannot be 
approached indefinitely for numerical reasons]^] Both posi- 
tive and negative travel-time shifts are seen. The effect is 
quite small for deeper waves (smaller i) , where the wave ap- 
proaches the surface almost vertically and hence interacts 
only weakly with vertical magnetic field, but is greatly en- 
hanced for shallow skippers (|zi| < 5 Mm). 

3.2.2 3D Inclined Magnetic Field 

A total of 3 x 4 x 3 = 36 cases are considered, corresponding 
to three magnetic field strengths B, 1 kG, 2 kG, and 3 kG; 
four lower turning depths zi, —3 Mm, —5 Mm, —10 Mm, 
and —15 Mm; and three wave frequencies, 3 mHz, 4 mHz, 
and 5 mHz. For z\ — — 3 Mm we adopt Zb = — 2 Mm and 
z s = —1 Mm; for z\ = —5 Mm we set Zh = —4 Mm and 
z s = —2 Mm; and for the other two cases Zb = —8 Mm and 
z s = —4 Mm. The top boundary conditions are applied at 
Zt = 2 Mm throughout. The spherical harmonic degree £ 
for each case is listed in Table [l] These are typical values 
in local helioseismology. Note that no attempt is made to 
place these cases on p-mode ridges, in recognition of the 
"scattering experiment" nature of our calculations: a wave 
is "fired" from Zb and examined when it reflects back to Zb. 
It is not necessary for it to be a normal mode, as the bottom 
boundary condition is not homogeneous]^] 

Figure [4] shows the fast wave reflection coefficient 1Z and 
phase shift Acp for B = 1 kG. The twelve constituent frames 

2 The governing differential equations for the magnetic case are 
either fourth order (2D) or sixth order (3D), whereas the non- 
magnetic equation is second order (see Appendix Ja|. The three 
cases therefore require different numerical solution processes, with 
different numbers of boundary conditions. As B — > 0, the slow 
and Alfven wavelengths vanish, so letting B become very small 
requires ever decreasing step-lengths, which of course cannot be 
continued through to the limit. 

3 Strictly speaking, there are no normal modes of the mag- 
netic solar model, as magnetic wave leakage produces pseudo- 
modes instead, resulting in complex eigenfrequencies or eigen- 
wavenumbers (ICally, Bogdan, & Zweibel||1994| iBogdan & Cally! 
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Figure 2. Two-way travel-time perturbations as a function of 
frequency for the MMMS model with magnetic field suppressed, 
but its thermal consequences retained. The top three curves cor- 
respond to B = 3 kG, the middle three to 2 kG, and the bottom 
three to 1 kG. The full curves are for waves with bottom turning 
points at Z\ = — 15 Mm, the dashed curves for z\ = —10 Mm, 
and the dotted curves for z\ = — 5 Mm. 
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B (kG) B (kG) 

Figure 3. Travel-time perturbation At = —u)~ 1 Aip (seconds) for 
vertical magnetic field and four lower turning depths z\ (labelled) 
as a function of magnetic field strength (kG). Full curves: 3 mHz; 
dashed curves: 4 mHz; dotted curves: 5 mHz. 



and longer path length naturally results in positive travel- 
time shifts. 

These results are not definitive. Different thermal mod- 
els can produce quite different timings. For comparison, an 
alternate thermal model with reduced surface density and 
negative non-magnetic travel-time perturbations is briefly 
addressed in Section [3. 3| However, the point here is to test 
whether thermal and magnetic effects are additive, rather 
than to produce a fully realistic sunspot model. 



3.2 Magnetic Anomaly Only 

The effects of a uniform magnetic field superimposed on un- 
modified Model S are examined next. These are presented 
in detail, as a major theme of this paper is to emphasize the 
direct magnetic effects which are too-often ignored in local 
helioseismology. 



|1997| |Crouch fc Cally||2005[ l. However, the sub-surface acoustic 
part may be decoupled from the magnetic in a perturbation sense 
Cally |20 05 ) , or by using the divergence of the displacement as is 
done here. 
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Figure 4. Fast wave reflection coefficeint "R, (vertical axis) and relative phase jump Ay? (colour shading) as a function of magnetic 
field inclination to the vertical 9 and orientation from the x-z plane </>. The twelve cases depicted are the combinations of four different 
turning depths Z\ (—3 Mm, —5 Mm, —10 Mm, and —15 Mm) and three different frequencies (3 mHz, 4 mHz, and 5 mHz), as labelled. 
The magnetic field strength is 1 kG throughout. The remainder of the flux is variously lost to the overlying atmosphere as acoustic or 
magnetic waves, or downward as magnetic waves. The colour shading represents Ay? on a periodic colour scale (see the legend disk at 
the top), with red corresponding to Ay? = 0° = ±360°, yellow-green to 90°, etc. A 4° X 10° 9-(f> grid was used throughout. 



correspond to the four turning depths and three frequencies, 
as labelled, and cover -84° < 6 < 84°, < <j> < 90°. Of 
course, results are independent of cj> at 9 = 0°, and Ay? is 
small there in agreement with Figure [3] However, substan- 
tial phase perturbations are found at large inclination 9, 
especially at 5 mHz. Clearly, fast mode reflection TZ is most 
pronounced at large 9, corresponding to penumbra. This is 
to be expected, as mode conversion is strongest when the 
wave strikes the magnetic field at z cq at large attack angle 
(the angle between the wavevector and the magnetic field 



lines; see Schunker fc Cally|2006 l. At low frequency though, 
there is still substantial fast wave reflection throughout. 

As might be expected, Figures [5] and [6] show a gener- 
ally increasing phase perturbation with increased magnetic 
field strength. In most cases, maximum phase perturbation 
is attained at high inclination 9, which is also where the re- 
turning fast wave flux is highest. (Of course, not all parts of 
all frames are relevant to the Sun. Highly inclined field for 
instance represents penumbra, where field strengths will be 
somewhat less than 2 kG generally.) 

Figure [7] illustrates that travel-time perturbation At 
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Figure 5. As for Figure]?] but with 2 kG magnetic field. 



at large field inclination varies essentially linearly with in- 
creasing B, and also increases in magnitude with increasing 
frequency. The phase perturbation is somewhat reduced at 
= 90° compared with = 0°. 

It is striking that phase perturbation at high field in- 
clination is relatively insensitive to turning depth z\. It is 
also notable that the travel-time perturbation is uniformly 
negative here. Both features are in stark contrast to the ver- 
tical field case. The probable cause is that the attack angle 
is large for all cases in the near-horizontal field, whereas it 
becomes quite fine for deeper rays in vertical field. 

A prominent characteristic of Figures[4}|6j especially the 
last two, is the banded structure of the phase perturbation 
shading, and the associated "W-shaped" return flux vari- 
ation with 9. There is a sharp distinction between Aip at 
small and large inclination, with low fast mode reflectivity 



1Z typically found at the transition region between the two. 
We focus on the flux first. 

The complex outgoing flux structure with varying field 
inclination 9 of one of the 5 mHz 3D cases is illustrated 
m Figure [8] The full curve shows the downgoing acous- 
tic flux at 2b, clearly displaying very distinct "central" 
and "wing" regions. The reason for this is clear. The up- 
ward acoustic flux (dashed curve) turns on sharply once 
\9\ > arccos(o)/tj c ) = 16.2° in this case. The upward acous- 
tic losses are maximal where the attack angle is smallest, 
around 30° here, in agreement with the complex eigenmode 
calculations of |Crouch fc Cally| p005] ) . They are small for 
negative 9 because the attack angle is large there. On the 
other hand, the downward magnetic losses (dotted curve) 
display the opposite behaviour. This is because they result 
from mode transmission on the downward leg of the fast 
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3 mHz, z, = -3 Mm, B = 3 kG 




Figure 6. As for Figure]?] but with 3 kG magnetic field. 



ray, where now attack angle is smallest around —30° and 
large at positive 8. The upgoing Alfven flux (dot-dashed 
curve) also displays an inclination asymmetry. So, in this 
and many other cases, the "centre/wing" dichotomy in 1Z is 
easily understood from flux graphs, and is in complete ac- 
cord with the generalized ray-theoretic insights of Schunkcr 



& Cally (20061. This angle-dependent "absorption" of helio- 



seismic waves by sunspots also explains the Hankel analysis 
results of 



Braun 



( 1995 1 (see also Cally, Crouch, & Braun 



theoretic description of phase jumps at mode conversion 
points. However, the behaviour is very reminiscent of the 
"central dips" seen in phase jumps across mode conversion 



levels in Figure 8 of Cally (20081. We therefore tentatively 



attribute the effect to mode-conversion rather than turn- 
ing point phase jumps, which the isothermal model suggests 
should be relatively insensitive to field inclination. 



2003| |Crouch et al.|2005| where calculated complex eigenval- 3.3 Combined Magnetic and Thermal Anomalies 

ues are used to reproduce the Hankel absorption and phase T . , , . , „ , . 

, . . , \ mi , , n r il i Let us now explore the combined effects of the magnetic 

shift data quite accurately), lhe weak dependence of ab- . ,. . ,. „. . ,. , ,. , 

_ . „ v _ , . . . and thermal anomalies, lhe travel-time perturbation due to 

sorption coefficient a = 1 — TZ on d> was previously noted . . . . . .. . , 

both magnetic and thermal anomalies may be represented 

by At(B, AT). Is this well approximated by the sum of the 

The reason for the banded phase structure is less ob- individual effects At(B,0) + Ar(0, AT)? Figure [9] makes it 

vious, in the current absence of a simple generalized ray- clear that the answer is "yes" for weak magnetic field (1 



by Crouch fc Cally ( 2005 1 
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Figure 7. Travel-time perturbations At = — 1 A(/p as a function of magnetic field strength (kG) at high field inclination, 9 = 84°, for 
turning depths z\ = —5 Mm (left column); z\ = — 10 Mm (centre column); and z\ = — 15 Mm (right column). The top row corresponds 
to the 2D case <j> = 0°, and the bottom row to the perpendicular case <j> = 90°. The full curve corresponds to a frequency of 3 mHz 
throughout, dashed to 4 mHz, and dotted to 5 mHz. Highly inclined field in sunspot pcnumbrae would not normally exceed about 1.5 
kG, so travel-time advances of around 1 minute may be expected from these results. 
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Figure 9. Travel-time perturbations as a function of magnetic field angle 9 with = 0° for frequencies of 3 mHz (top row) and 5 mHz 
(bottom row) and 1 kG magnetic field (left column), 2 kG (centre column) and 3 kG (right column). The full curve represents At(B, AT), 
the long-dashed curve Ar(B,0) + At(0, AT), short-dashed Ar(B,0), and dotted At(0, AT). The At range plotted corresponds to 
— 180° ^ Aip ^ 180°, and in several cases displays a 360° "wrap-around". 



kG, left column), where the full and long-dashed curves are 
very close, but a resounding "no" at 3 kG, unless the field 
is nearly vertical. At 2 kG, we have a qualified "yes" for 
6*| < 25°. 

Of course though, highly inclined field (penumbra) typ- 
ically does not exceed 2 kG, so the large curves at 
B — 3 kG are moot. We incorporate the tendency for 
field strength to decrease with field inclination by setting 
B = B max exp[-(r?/90°) 2 ], with B max = 3 kG, which yields 



1.1 kG at 6 = 90°. Figure [To] then displays rather com- 
plex behaviour. Negative travel-time shifts appear charac- 
teristic of the penumbra, though positive shifts are typi- 
cally seen in the umbra (small |#|), in accord with the ob- 
servations of Couvidat & Rajaguru (20071. In the <j) = 0° 



cases (2D), this is due to the 360° ambiguity in Aip, both 
for deep (left column) and shallow (right colum) skipping 
waves. For the transverse 3D case {<j> = 90°, centre col- 
umn), all solutions displayed stay within the principal win- 
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Figure 8. Outgoing wave energy fluxes for the 5 mHz case with 
z\ = — 10 Mm, B = 2 kG, </> = 30°, corresponding to a slice across 
the last graph in the third row of Figure [5] (though with improved 
1° resolution). Full curve: F a b = TZ. (outward acoustic flux at the 
bottom - the "helioseismic" field); Dashed curve: F a t (outward 
acoustic flux at the top); Dot-dashed curve: F m t (Alfvcnic flux 
at the top); and Dotted curve: F m i, (combined slow and Alfvcn 
magnetic flux at the bottom). Of course, F a b + F a t + F m t + F m j, = 
Fi n = lj where F ln is the injected acoustic flux. 



dow. For the most part, the assumption of additive ther- 
mal and magnetic effects is rather poor in the penumbra, 
and only sometimes good in the umbra. In fact, it appears 
that At(B, AT) rs At(B, 0) is a better approximation than 
At(B, AT) w At(B, 0) + Ar(0, AT) in penumbrae. 

Finally, the calculations of Figure [10] are repeated with 
an alternate heuristic thermal model in which the density 
is not held fixed. In fact, the scale height H decreases in 
the near-surface layers as they are cooled, sufficiently fast 
that w c = (c/2H)\/l — 2H' actually increases. This results 
in shorter travel paths and negative travel time shifts in the 
non-magnetic case, despite the reduced sound speed. Fig- 
ure [TT] illustrates the results. Again, it is clear that direct 
magnetic effects dominate At in penumbrae, but the um- 
bra! shifts are predominantly thermal. However, we might 
also conclude that At(B, AT) w At(B, 0) + At(0, AT) does 
rather better in this model, probably because the magnitude 
of Ar(0, AT) is typically much smaller. 



4 DISCUSSION 



The analysis in Cally ( 2008 ) largely proceeds by following a 
fast mode ray through the surface layers as it first propagates 
upward, refracts from the Alfven speed gradient, and finally 
propagates downward. The ray is used as a "sampler" of the 
known exact wave field in the 2D isothermal case: i.e., de- 
termining the exact phase at the ray-theoretic position x(t) 
at the time t = T when the ray returns to z = z\,. One may 
also compare this to the ray-theoretic phase S(T) calcu- 
lated along with the ray path. We have chosen not to apply 
that procedure here, largely because the ray integration is 
unreliable in realistic solar models near the surface, for rea- 
sons associated with the small vertical length scales, and the 
definition of the acoustic cutoff frequency (see |Schunker fc| 
Cally 2006). On the other hand, the definition of Sip in Equa- 
tion (|2j) samples phase "in place", i.e., 5p = — arg(C-/C+) 
rather than — arg(C_/C+) + kx(T) — uoT. It pays no heed 
to the differing ray paths and travel times. Nevertheless, 



the "in place" Sip is made useful by comparison to the non- 
magnetic phase difference between downgoing and upgoing 
acoustic waves at Zb, i.e., Aip — Sip — Spo. This is consis- 
tent with the standard local helioseismic perspective, where 
comparison is always made to quiet Sun values. 

A clear lesson to be drawn from all the magnetic-only 
cases presented in Section |3.2| is that phase perturbation 
Aip is most pronounced for highly inclined magnetic field, 
characteristic of sunspot penumbrae. The significance of this 
correspondence is enhanced by the maximal fast mode re- 
flection coefficient 1Z at large inclination. Therefore, penum- 
bra (and perhaps plage with its predominantly horizontal 
overlying canopy field) may be expected to be most reflec- 
tive and also most productive of phase astigmatism. This 
appears to be consistent with observational evidence from 
TP and HH (|Couvidat fc Rajaguru|2007| pndsey fc Braun 



200r5}|Schunker et al.|2005| 

The travel-time perturbations presented here, in Figure 

3 (umbra) and Figure [7] (penumbra), and indirectly Figures 

4 [6] should be understood as due to direct magnetic effects 
on the waves, through mode conversion, transmission, and 
fast wave reflection. 

In addition to these effects though, the magnetic field of 
sunspots significantly alters the near-surface thermal struc- 
ture, and therefore the sound speed. This is addressed in Sec- 
tionpOwith an ad hoc thermal model (MMMS). We find a 



definite tendency for negative travel-time shifts of around a 
minute in highly inclined (and comparatively weak) penum- 
bral field. The behaviour in more vertical (and stronger) 
umbral field is rather complex, often associated with the 
360° phase ambiguity, but is typically positive if the window 
— 180° ^ Aip ^ 180° is assumed. This raises a further doubt 
about the wisdom of interpreting At — —cj~A(p as a true 
"travel-time" perturbation: the 360° wrapping can turn an 
actual large negative travel-time shift into an apparent pos- 
itive perturbation quite easily. It is tempting to hypothesize 
that this may be responsible for the positive umbral shifts re- 



ported by Couvidat & Rajaguru (20071. Alternate plausible 



thermal models may produce quite different thermal shifts 



(figure 111, but the general conclusions are unchanged. 



Based on the results illustrated in Figures [T0| and [TT] a 
crude rule-of-thumb appears to be that travel-time pertur- 
bations in umbrae are predominantly thermal, whereas in 
penumbrae they are mostly magnetic. This is a neat explana- 



tion for the umbra/penumbra dichotomy identified by Cou- 



vidat & Rajaguru (20071. The explanation for this behaviour 



is clear: for heliosesimic rays propagating steeply to the sur- 
face, the direct interaction with the near-vertical magnetic 
field is small, resulting in substantial mode transmission 
(acoustic-to-acoustic) across the equipartition layer, almost 
as if the field were not there. On the other hand, mode con- 
version (acoustic-to-magnetic) dominates for highly inclined 
penumbral field. The resulting atmospheric fast wave quickly 
reflects back downward to rejoin the helioseismic wave field, 
but with a very different phase (and timing) to the non- 
magnetic case. 

Current local helioseismic techniques are particularly 
adept at measuring temperature, as well as flow speeds. 
However, they are easily confused by the intrinsic mag- 
netic perturbations. It is hoped that the results presented 
here will point the way to interpreting observed local- 
helioseismic data in sunspots. Further forward modelling 
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Figure 10. Travel-time perturbations in the MMMS model as a function of magnetic field angle 8 with z\ = — 15 Mm and <j> = 0° (left 
column), z\ = —15 Mm and <j> = 90° (centre column), and Z\ = —5 Mm and <j> = 0° (right column), for B = 3 exp[— (0/9O°) 2 ] kG 
and three frequencies: 3 mHz (top), 4 mHz (centre), and 5 mHz (bottom). The full curve represents At(B, AT), the long-dashed curve 
At(B, 0) + At(0, AT), short-dashed At(B, 0), and dotted Ar(0, AT). The At range plotted corresponds to -180° ^ Aip < 180°, and 
in several cases exhibits a 360° "wrap-around" . 



based on more or less realistic magnetic and thermal sunspot 



models (Cameron, Gizon, fc Duvall||2008[ Moradi & Cally 



2008 Moradi, Hanasogc, & Cally 2008|) is also warranted. 
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Figure 11. Same as Figure [To| but for a thermal model with reduced density as well as temperature, thereby producing negative 
non-magnetic travel-time perturbations (dotted curves). 



Moradi, H., Cally, P.S., 2008, Solar Phys., 251, 309 
Moradi, H., Hanasoge, S.M., Cally, P.S., 2008, ApJL, ac- 
cepted 

Schunker, H., Braun, D.C., Cally, P.S., Lindsey, C: 2005, 

ApJL, 621, L149 
Schunker, H., Cally, P.S., 2006, MNRAS, 372, 55 
Weinberg, S., 1962, Phys. Rev., 126, 1899 



APPENDIX A: 
SOLUTION 



NONMAGNETIC AND WKB 



In the nonmagnetic case, the wave equation may be written 
in the form 



*" =u z Q(z)V, 
where 

uj 2 Q(z) = - 



I ^ -n 



(Al) 
(A2) 



and N is the Brunt-Vaisala frequency. The acoustic cutoff 
frequency is specified by lo 2 = (c 2 /4H 2 )(1 — 2H'), where H 
is the density scale height. (For numerical reasons associated 
with difficulties in calculating H' smoothly and reliably in 
tabulated solar models, this is not the formulation we use for 
numerical solution in the nonmagnetic case, though math- 
ematically it is equivalent to it.) Equation ( |A1[ ) is in con- 
venient form for WKB solution, with lo^ 1 taking the role 



of the perturbation parameter. The WKB solution in the 
acoustic cavity, where Q is negative, is 



exp 



(A3) 



as uj — > oo. See Cally (2008) for further details. A linear 
combination of /+ and /_ may be used to fit the numerical 
magnetic solution ^ at sufficient depth that the magnetic 
and acoustic waves have decoupled. This identifies the up- 
going and downgoing components. 



